
function PlotPseudoArc(yy , aux)
%
% 伪弧长延拓画图
%
% 作者：张晨
% 邮箱：chenzhang@csu.ac.cn
% 单位：中国科学院空间应用工程与技术中心，空间探索室
% 时间：2022年02月24日
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 画图
var1 = yy(1 : 4);
varf = yy(4 * (aux.nNode - 1) + 1 : 4 * aux.nNode);
tau1 = yy(4 * aux.nNode + 1);
tauf = yy(4 * aux.nNode + 2);

% 【新增】为提高收敛性，第一个节点（LEO）和最后一个节点（LLO）使用轨道根数表示
xx1 = Bcr4bp_var2xxEMRot2D(var1 , 'leo' , aux)';
xxf = Bcr4bp_var2xxEMRot2D(varf , 'llo' , aux)';

% 收敛轨道画图
options = odeset('Reltol' , aux.tol , 'AbsTol' , aux.tol);
[seg_tt , seg_xx] = ode113(@Bcr4bp_EqmEMRot2D , [tau1 , tauf] , [xx1] , options , aux);

% ---------------------------------- EMRot画图 ----------------------------------
figure(1);

% 画转移轨道
plot(seg_xx(: , 1) , seg_xx(: , 2) , 'LineWidth' , 1 , 'color' , [0.8 , 0.8 , 0.8]);

% ---------------------------------- SB1Rot画图 ----------------------------------
figure(2);
[tau_SB1Rot , xb2c_SB1Rot] = Bcr4bp_EMRot2SB1RotMtx2D(seg_tt , seg_xx , aux);

% 画转移轨道
plot(xb2c_SB1Rot(: , 1) , xb2c_SB1Rot(: , 2) , 'LineWidth' , 1 , 'color' , [0.8 , 0.8 , 0.8]);

end
